Discovery of a ZZ Ceti in the Kepler Mission Field 
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ABSTRACT 

We report the discovery of the first identified pulsating DA white dwarf in 
the field of the Kepler mission, WDJ1916+3938 (Kepler ID 4552982). This ZZ 
Ceti star was first identified through ground-based, time-series photometry, and 
follow-up spectroscopy confirm it is a hydrogen-atmosphere white dwarf with T e g- 
= 11,129 ± 115 K and log g = 8.34 ± 0.06, placing it within the empirical ZZ Ceti 
instability strip. The object shows up to 0.5% amplitude variability at several 
periods between 800 - 1450 s. Extended Kepler observations of WD J1916+3938 
could yield the best lightcurve, to-date, of any pulsating white dwarf, allowing 
us to directly study the interior of an evolved object representative of the fate of 
the majority of stars in our Galaxy. 

Subject headings: stars: white dwarfs-stars: individual (WD J191643. 83+393849. 7)- 
stars: oscillations (including pulsations)-stars: variables: general 



1. Introduction 



Hydrogen atmosphere (D A) white dwarf stars a ccount for more than 80% of all spec- 
troscopically identified WDs (jEisenstein et al.l 120061 ) . When these objects cool to around 
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12,500 K, the hydrogen in their non-degenerate atmosphere forms a partial ionization zone, 
impeding energy transport and leading to global g-mode oscillations. It is believed that all 
DA white dwarf stars wit h carbon/oxygen cores un dergo such pulsations when they reach the 
appropriate temperature (jCastanheira et al.ll2007l ). allowing us to probe a stellar population 
representative of the fate of most stars in our Galaxy, including our Sun. 

Asteroseismology provides a unique opportunity t o investigate the internal structure 
of these stars (see reviews by IWinget fc Keplerl 120081 ; iFontaine fc Brassard! 120081 ). Given 
enough observed modes, we can put accurate estimates on the overall stellar mass, hydrogen 
and helium layer masses, internal chemical transition zones, magnetic field strength, and 
rotation rate. Performing asteroseismology on a WD benefits from its degenerate interior, 
which makes its structure less complicated than for stars still undergoing fusion. 

In an effort to reduce aliasing caused by gaps in the data, much effort has been expended 
to obtain uninterrupted time-series observations of these v ariable DA (DAV) stars (also 
known as ZZ Ceti stars), using the Whole Earth Telescope ( iNather et al.lll990l ). However, 
these ground-based campaigns rarely extend beyond a few weeks, fundamentally limited by 
weather and moonlight, not to mention the availability of telescope time. 

Space-based photometry has revolutionized the extended coverage of stars. The Kepler 
spacecraft was launched in March 2009, and while its primary mission is to detect Earth- 
sized planets around Sun-like stars, its extensive lightcurves are transforming many fields 
in asteroseismology. The spacecraft stares at some 156,000 stars and has the capability to 
monitor 512 objects with one-minute cadence. Month-long observations typically yield a 
duty cycle of more than 97% with micromagnitude precision, depending on the brightness 
of the target. 



Despite an expansive search (lOstensen et al.ll2010l . l2011al ) in which 17 compact objects 
were targeted by Kepler as potential pulsating white dwarf stars, none were observed to 
vary above the 4cr detection limit. One DAV candidate, KIC10420021, was monitored for 7 
months and is stable to at least 60 parts-per-million (ppm). Kepler has monitored a uniquely 
variable DA, but t hat system is best expla ined by a magnetic spin-modulated WD model, 
not by pulsations (IHolberg &: Howell! 1201 ll ). However, as part of a recent auxiliary search, 
0stensen et~aD toil\h announced the successful detection of the first pulsating WD in the 
Kepler field, a V777 Her star. It was quite a fortuitous discovery, as this was just the second 
DB found in the field and has a Kepler magnitude of Kp = 18.46. 

Using facilities at the McDonald Observatory in west Texas, we have discovered the first 
ZZ Ceti star in the Kepler field, WD J191643.83+393849.7 (Kepler ID 4552982, hereafte r 
WD J1916+3938). This object has a Kepler magnitude of Kp = 17.85 fjBrown et al.ll201ll ). 
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Previous Kepler observations of such faint, blue objects have reached better than 130 ppm 
(0.013%) precision with just one month's worth of data. Extended short-cadence Kepler 
observations will provide a detailed look at the pulsation spectrum of this cool ZZ Ceti star. 



2. High-Speed Photometric Observations 
2.1. Selection and Time-Series Photometry 



Potential pulsating WDs in the Kepler field were selected f rom a recent list of hig h 
reduced proper motion stars in the SuperCOSMOS Sky Survey (IRowell &: Hamblyl 1201 ll ). 
From their catalog of roughly 10,000 candidate WDs, we found 20 that actually fell on 
Kepler CCDs. Five of those had already been surveyed by Kepler and found not to vary 
f l0stensen etaJboid . lioilah - 



One of the remaining 15, WD J1916+3938, fell in the footprint of the POSS-II survey, 
and showed (B—R), (R—I) colors extremely close to values for a typical ZZ Ceti star, so it 
was targeted for time-series photometric observations. The discovery observations were taken 
at the McDonald Observatory using the Argos instr ument, a frame-transfer CC D mounted 
at the prime focus of the 2.1m Otto Struve telescope (jNather &: Mukadamll2004j ). Data were 
obtained through a 1mm BG40 filter to reduce sky noise. 

Photometric observations were carried out over 6 nights in May 2011, for a total of more 
than 21 hours of coverage. Exposure times ranged from 5 to 15 seconds, depending on seeing 
and cloud conditions. 

We performed weighted, circular aperture photometry on the calibrated frames using the 
external IRAF package ccdJisp written by Antonio Kanaan. We divided the sky-subtracted 
lightcurves using at least 3 brighter comparison stars in the field to allow for fluctuat ions in 
seeing and cloud cover. Using the WQED software suite ([Thompson &: Mullallyll2009l ). we fit 
a second-order polynomial to the data to remove the long-term trend caused by atmospheric 
extinction, and applied a timing correction to each observation to account for the motion of 
the Earth around the barycenter of the solar system. 

Figure [1] shows a representative lightcurve, smoothed by a four-point moving average 
window for visualization purposes with the brightest comparison star offset below. 
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2.2. Ground-Based Lightcurve Analysis 

We performed a Fourier transform of the entire data set, 7,454 points spread over nearly 
a month for a duty cycle of 3.1% (Figure [2]). This FT extends to the Nyquist frequency 
for Kepler short-cadence observations; there is nothing above the 4<j threshold for higher 
frequencies. The window function is quite messy, as there are variably spaced gaps in the 
data. 

While Kepler will allow us to improve our determination of the periodicities present 
in this star, we attempt to identify the periods present in all our May 2011 data. Tabled] 
lists the highest seven peaks, all with power above four times the mean FT level. The 
large aliasing present makes it nearly impossible to determine, with much certainty, the true 
periods in this short, ground-based data set, so we have omitted the formal least-squares 
errors. It appears this star has several excited modes between 800 - 1450 s, although the 
amplitudes appear unstable from night to night. There is consistently power well above the 
noise around 970 s and 820 s each night, although the level of excitation appears inconsistent. 

For comparison sake, we also took an FT of all the data for the brightest comparison 
star in the field. The average amplitude over all frequency space is 0.03%, the same level as 
our formal least-squares amplitude errors. 

3. Spectroscopic Analysis 

3.1. Spectroscopic Observations 

Excitement among the Kepler compact object community, specifically Working Group 11 
within the Kepler Asteroseismic Science Consortium (KASC), led to multiple spectroscopic 
follow-ups, and WD J1916+3938 was observed using 4 different telescopes, especially after 
our first spectrum from the Nordic Optical Telescope (NOT) confirmed the object as a DA 
within the ZZ Ceti instability strip. Fits to these observations were used to derive an effective 
temperature and surface gravity. 

Observations using the 9.2m Hobby-Eberly Telescope (HET) at McDonald Observatory 
were carried out with the Marcario Low Resolution Spectrograph (LRS) using its g2_2.0 
setup. The observations using the 2.6m NOT at Roque de los Muchachos Observatory were 
done in low-resolution with the ALFOSC spectrograph and its grism #7. The observations 
using the 3.5m telescope at Calar Alto Observatory (CAHA) were done in medium-resolution 
with the TWIN spectrograph, using grating T06 for the red arm and grating T12 for the 
blue arm. Finally, the observations using the 4m Mayall Telescope at Kitt Peak National 
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Observatory (KPNO) were done in low-resolution with the R-C CCD Spectrograph. Each 
spectral image and flux standard were reduced using standard long-slit IRAF routines. 

The KPNO spectrum of WD J1916+3938, showing the Balmer lines FL3 through H 9 , can 
be seen in Figure |3j Overplotted is the best-fit DA model spectrum. 

We have fit model grid s pectra to each of these observations, using the procedure out- 
lined in lBergeron et al.l ( 119921 ). The WD models used for the spe ctroscopic fitting were kindly 
provided by Detlev Koester and are described in iKoesteri J201(ih_ Balmer lines in the models 
were calculated with the modified Stark broadening profiles of iTremblay &: Bergeron! ( 120091 ). 
kindly made available by the authors, and use the mixing length prescription ML2/a=0.6. 
Table [2] lists the best fits to the observations. Using a weighted mean of all four determina- 
tions, we adopt a best value of T e g = 11,129 ± 115 K and log g = 8.34 ± 0.06. This places 
the object within the empirical ZZ Ceti instability strip. 



3.2. Implications 



ZZ Ceti stars at the hotter, blue edge of the instability strip tend to pulsate with short 
period (100 - 300 s) modes that have shown to be extremely stable in period and amplitude. 
Such stars have be en the focus of wo r k to detect an evoluti onary rate of period change of 
these modes (e.g., Kepler et al.l 120051 ; iMukadam et al.l 120091 ). However, the ZZ Ceti stars 
near the cooler, red edge of the instability strip with effective temperatures near 11,000 
K typically feature much longer periods (> 600 s) and have more excited modes. These 
pulsations are typically unstable in amplitude and phase. 

Spectroscopic fits show that WD J1916+3938 has an effective temperature that places 
it in the cooler region of the ZZ Ceti instability strip. This result is duly confirmed by our 
time-series photometric observations. WD J19 16+3938 shows long-period oscillations — its 
weighted mean period (IMukadam et al.l 120061 ) lies around 1,003 s — that do not appear 
stable in amplitude. For example, during the runs from May 27-29, the power observed 
in the mode near 824 s went from (0.66 ± 0.06)% to (0.88 ± 0.05)% to (0.43 ± 0.05)%, 
respectively. 

Many modes in other cool ZZ Ceti stars have amplitudes that grow and decay in a way 
that cannot be explained by beating (e.g., iKleinman et al.lll998l ). so we expect the same 
is true with this star. The opportunity for nearly continuous Kepler observations of this 
amplitude modulation may afford us a completely new insight into the growth and decay 
of mode amplitudes, opening a new window into how these objects behave on a variety of 
timescales inaccessible from the ground. 
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Discussion and Conclusion 



We have discovered just the second pulsating white dwarf in the Kepler field, and the 
first ZZ Ceti star suitable for Kepler monitoring. Several modes appear to be excited from 
800 - 1450 s, with a fractional amplitude of up to 0.5%. Spectroscopic fits yield an effective 
temperature of 11,129 ± 115 K and a surface gravity of log g = 8.34 ± 0.06, which corre- 



sponds to a mass of 0.82 ± 0.04 M©from the evolutionary models of IWoodl ( 1l990l ). This 
spectroscopic temperature puts it on the cool, red edge of the ZZ Ceti instability strip, where 
pulsations may begin shutting down. 

Submitted for Kepler observations through KASC, the spacecraft is scheduled to com- 
mence short-cadence observations on this ZZ Ceti pulsator in October 2011 (Qll). By virtue 
of its impressive duty cycle, we anticipate that extended Kepler observations will yield the 
best lightcurve, to date, of any pulsating white dwarf, even those that are more than four 
magnitudes brighter. 

The lightcurves can be analyzed from a variety of asteroseismic angles. With enough 
modes, we will be able to put an accurate constraint on the overall mass of the WD. Per- 
haps more importantly, we may also extract an accurate estimate for the hydrogen layer 
mass, which controls the rate at which this star cools. This hydrogen mass fraction has b een 
constra ined through asteroseismolog y in fewer than a dozen DA WDs (IDolez et al.ll2006l ). al- 
though |Castajih^ira^jie^lerJ (120091 ) attempted to derive the parameter from the few modes 
present in an ensemble of 83 ZZ Ceti stars. Uncertainties in this parameter is imprinted as 
scatter in the white dwarf luminosity fu nction, which is used in cosmo chronolo gy to put an 
estimate on the age of the Galaxy (e.g., 



Winget et al.lll987t lHarris et al.ll200q ). 



The sensitivity of Kepler observations may also allow us to see a wealth of low-amplitude 
multiplet structure. Previous ground-based campaigns on cool ZZ Ceti stars, such as G29 - 
38, have revealed instability in the amplitude of these multiplets (IKleinman et al.lll998l ). 
Extended Kepler coverage will allow us to observe how the amplitudes of all modes in this 
star change with time. 

ZZ Ceti stars often exhibit so me nonlinear combination frequencies, which are believed to 
arise from convective driving (e.g. . I Br i ckhi 111 1 1 9 9 2| ; IWull200l[ ) . If observable, these combination 



frequ encies may be used to identify the spherical degree of the parent modes ( jYeates et al. 



20051 ). With an ideal K epler lightcurve, we would be able to compare this mode identification 
with the £ values as determined from the period spacings. Pulse-shape fitting of properly 
identified modes may also be used to constrain some of the convective p arameters in this 
star, especially the depth of the convection zone ([Montgomery et al.ll2010l ). 



For now, though, we eagerly await the unrivaled staring competition Kepler will have 
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with this cool ZZ Ceti star. 

We thank KASC WG11 for their enthusiasm, especially regarding spectroscopic follow- 
ups. This work is supported by the Norman Hackerman Advanced Research Program, under 
grants 003658-0255-2007 and 003658-0252-2009, by a grant from the NASA Origins Pro- 
gram, NAG5-13094, and by the National Science Foundation, under grant AST-0909107. 
The research leading to these results has received funding from the European Research 
Council under the European Community's Seventh Framework Programme (FP7/2007- 
2013) /ERC grant agreement n°227224 (PROSPERITY), as well as from the Research Coun- 
cil of K.U.Leuven grant agreement GOA/2008/04. This work was partly based on observa- 
tions made with the Nordic Optical Telescope, operated on the island of La Palma jointly 
by Denmark, Finland, Iceland, Norway, and Sweden. We thank the McDonald Observatory 
staff for their support, especially Dave Doss and John Kuehne. 

REFERENCES 

Bergeron, P., Saffer, R. A., & Liebert, J. 1992, ApJ, 394, 228 

Bergeron, P., Fontaine, C, Billeres, M., Boudreault, S., & Green, E. M. 2004, ApJ, 600, 404 
Brickhill, A. J. 1992, MNRAS, 259, 519 

Brown, T. M., Latham, D. W., Everett, M. E., & Esquerdo, G. A. 2011, arXiv:1102.0"342l 

Castanheira, B. G., et al. 2007, A&A, 462, 989 

Castanheira, B. G., & Kepler, S. O. 2009, MNRAS, 396, 1709 

Dolez, N., et al. 2006, A&A, 446, 237 

Eisenstein, D. J., et al. 2006, ApJS, 167, 40 

Fontaine, G., & Brassard, P. 2008, PASP, 120, 1043 

Harris, H. C, et al. 2006, AJ, 131, 571 

Holberg, J. B., & Howell, S. B. 2011, AJ, 142, 62 

Kawaler, S. D. 1987, Stellar Pulsation, 274, 367 

Kepler, S. O., et al. 2005, ApJ, 634, 1311 



-8 - 

Kleinman, S. J., et al. 1998, ApJ, 495, 424 

Koester, D. 2010, Mem. Soc. Astron. Italiana, 81, 921 

Montgomery, M. H., et al. 2010, ApJ, 716, 84 

Mukadam, A. S., Montgomery, M. H., Winget, D. E., Kepler, S. O., & Clemens, J. C. 2006, 
ApJ, 640, 956 

Mukadam, A. S., et al. 2009, Journal of Physics Conference Series, 172, 012074 

Nather, R. E., Winget, D. E., Clemens, J. C, Hansen, C. J., & Hine, B. P. 1990, ApJ, 361, 
309 

Nather, R. E. & Mukadam, A. S. 2004, ApJ, 605, 846 
Cstensen, R. H., et al. 2010, MNRAS, 409, 1470 
Cstensen, R. H., et al. 2011a, MNRAS, 414, 2860 
Cstensen, R. H., et al. 2011b, ApJ, 736, L39 

Rowell, N. & Hambly, N. C. 2011, arXiv:astro-ph.GA/1102.3193v2 

Thompson, S. E., & Mullally, F. 2009, Journal of Physics Conference Series, 172, 012081 

Tremblay, P.-E., & Bergeron, P. 2009, ApJ, 696, 1755 

Winget, D. E., Hansen, C. J., Liebert, J., van Horn, H. M., Fontaine, C, Nather, R. E., 
Kepler, S. O., & Lamb, D. Q. 1987, ApJ, 315, L77 

Winget, D. E., & Kepler, S. O. 2008, ARA&A, 46, 157 

Wood M. A., 1990, PhD thesis, Univ. Texas , Austin 

Wu, Y. 2001, MNRAS, 323, 248 

Yeates, C. M., Clemens, J. C, Thompson, S. E., & Mullally, F. 2005, ApJ, 635, 1239 



This preprint was prepared with the A AS IATgX macros v5.2. 



-4 



Comparison 
J L_ , , I , , , L 



2000 4000 6000 

Time (s) 



Fig. 1 — A portion of the lightcurve for WD J1916+3938, taken on 2011 May 11, that has 
been smoothed by a four-point moving average. The red, offset lightcurve is of the brightest 
comparison star in the field. 
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Fig. 2. — A Fourier transform of our May 2011 data on WD J1916+3938. The red transform 
has been computed after prewhitening by the seven periodicities listed in Table dj which 
removes all peaks above four times the mean FT level, a. We have ignored the low-frequency 
peaks in this FT, with periods longer than 2,000 s, which are likely to be noise from variability 
in the Earth's atmosphere. 
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Fig. 3. — The KPNO spectrum of WD J1916+3938, showing the Balmer lines Hg through 
H 9 . Overp lotted in red is the best-fit DA model spectrum. 
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Table 1: Frequencies present in WD J1916+3938 from May 3-31. 



ID 


Freq. (/xHz) 


Period (s) 


Amp. (%) 


h 


1032.1 


968.9 


0.44 


h 


1213.7 


823.9 


0.38 


h 


1070.1 


934.5 


0.36 


u 


1198.9 


834.1 


0.32 


h 


918.3 


1089.0 


0.25 


h 


696.0 


1436.7 


0.24 


h 


854.8 


1169.9 


0.23 



Table 2: Spectroscopically determined parameters. 



Instrument 


Date Observed 


Wavelength Range (A) 


T cS (K) 


log# (dex) 


NOT/ALFOSC 


2011-05-22 


3850-6850 


11,152 ± 80 


8.24 ± 0.06 


C AH A/TWIN 


2011-06-17 


3500-5450/5700-6790 


11,116 ± 139 


8.36 ± 0.09 


HET/LRS 


2011-06-14 


4300-7200 


11,104 ± 150 


9.51 ± 0.43 


KPNO/RCS 


2011-06-09 


3650-5120 


11,115 ± 92 


8.43 ± 0.05 


Adopted 






11,129 ± 115 


8.34 ± 0.06 



